Exploration of the Immuno-Inflammatory Potential Targets of Xinfeng Capsule in Patients with Ankylosing Spondylitis Based on Data Mining, Network Pharmacology, and Molecular Docking

Objective This study aimed to ascertain the immuno-inflammatory molecular targets of Xinfeng capsules (XFC) in the treatment of ankylosing spondylitis (AS) based on data mining, network pharmacology, and molecular docking. Methods The efficacy of XFC in the treatment of AS was assessed by clinical data mining. Network pharmacology was utilized to establish a network of the targets for XFC active ingredients in the treatment of AS. The binding mode and affinity of XFC active ingredients to the key targets for AS were predicted using molecular docking. Results XFC significantly diminished immuno-inflammatory indicators of AS. In total, 208 targets of XFC were obtained from the TCMSP database and 629 disease targets of AS were screened from the GeneCards database, which were intersected to yield 57 targets of XFC in the treatment of AS. Protein-protein interaction, gene ontology, and Kyoto genome encyclopedia analyses showed that XFC might activate TNF and NF-κB signaling pathways. Quercetin, kaempferol, triptolide, and formononetin had free binding energies < -9 kcal/mol to inflammatory targets (TNF and PTGS2) in the molecular docking analysis of XFC-active ingredients, indicating that TNF and PTGS2 might be the targets of the action of XFC. Conclusions Collectively, XFC had a significant therapeutic effect on AS. Specifically, the active ingredients of XFC, including quercetin, kaempferol, triptolide, and formononetin, inhibited the inflammatory response in AS by downregulating TNF and PTGS2 in the TNF and NF-κB signaling pathways.


Introduction
Ankylosing spondylitis (AS) is a kind of chronic inflammatory arthritis that mainly afflicts the axial and sacroiliac joints [1]. It occurs more frequently in men with a high incidence and disability rate [2,3]. If AS is not treated appropriately, this inflammatory damage normally causes chronic back pain and morning stiffness in patients, thereby resulting in the restriction of spinal mobility and functional impairment status of patients [4][5][6], which seriously affect the life quality of patients. According to existing research, the pathogenesis of AS is related to many factors. However, the specific etiology is still unclear [7].
At present, the most common clinical treatments for AS mainly include symptomatic treatments to relieve the pain of patients. However, the effects of these treatments are not ideal and often lead to certain side effects [8]. erefore, the development of new AS treatment strategies is essential.
Xinfeng capsules (XFC), as an in-hospital preparation of the First Affiliated Hospital of Anhui University of Traditional Chinese Medicine (TCM), has the effects of invigorating spleens, replenishing Qi, removing dampness and numbness, and activating blood and collaterals [9]. e previous study has proved the safe and controllable curative effects of XFC in the clinical application of AS [10].
Data mining technology is a technology that explores potentially valuable knowledge from a large amount of data, which has been widely used in the medical domain [11]. Retrospective data mining research can find meaningful results from clinical data. Network pharmacology technology has been utilized to predict the mechanism behind the action of drugs by constructing a complex network among drugs, ingredients, targets, diseases, and pathways, which is suitable for the studies of complex system mechanisms of TCM [12].
is piece of research sets out to screen the chemical active ingredients and disease-related targets of XFC in AS treatment on the basis of the methods of clinical data mining, network pharmacology, and molecular docking through data analysis under the guidance of bioinformatics theory to assess the impacts of the treatment process. e enrichment analysis of the highly correlated targets was conducted using molecular docking technology to fit the binding activity between the drug-active ingredients and target protein molecules to explore the potential effective ingredients and molecular mechanism of XFC in the treatment of AS and to further clinically verify the mechanism underlying the action of XFC.

Acquisition of Clinical Data and Analysis of Association
Rules.
e inpatient data were collected for 908 patients who were AS inpatients between July 2009 and December 2020 in the Department of Rheumatology and Immunology of the First Affiliated Hospital of Anhui University of Chinese Medicine. e data included the usage of XFC and the immuno-inflammatory indexes, such as the erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), immunoglobulin A (IgA), immunoglobulin M (IgM), immunoglobulin G (IgG), complement component 3 (C3), and C4.
is study was approved by the Medical Ethics Committee of Anhui Provincial Hospital of Traditional Chinese Medicine (2020AH-08). All study participants provided the written informed consent. e Apriori module in the IBM SPSS Modeler 18.0 software was adopted to analyze the correlation between XFC and immune-inflammatory indicators. e formula is as follows [13]: (1)

Acquisition and Screening of Pharmaceutical
Ingredients. e chemical active ingredients of Astragalus, coix seed, Tripterygium wilfordii, and centipede in the XFC prescription were acquired from the traditional Chinese medicine system pharmacological analysis platform (TCMSP, http:// tcmspw.com/) [14]. e above ingredients were retrieved as per the toxic pharmacokinetic parameters (absorption, distribution, metabolism, and excretion) [15] with the screening criteria of the oral availability of the drug ≥30% and the drug similarity ≥0.18 [16,17]. e relevant literature was employed to improve the data.

Prediction of Drug Component Targets and Acquisition of Disease-Related Targets.
e targets of the active ingredients of XFC were obtained using the predicted targets of the active ingredients in drugs from the TCMSP database. e predicted targets of the chemical active ingredients were harvested after the completion of the standardization conversion in the Universal Protein Resource (UniProt). en, the network visualization software, Cytoscape [18], was applied to construct and analyze the network characteristics and the targets of the active ingredients to clarify the relationship between the active ingredients and the targets of XFC.
With "AS" utilized as a search term, the disease-related targets were screened in the databases of GeneCards (https:// www.genecards.org/) and OMIM (http://www.omim.org/). e data of these two databases were intersected, and the unique values were retained to obtain known disease-related targets of AS.

Protein-Protein Interaction (PPI) Network Construction and Analysis.
e overlapping targets were imported into the STRING [19] protein interaction database (http:/stringdb.org/) with the species defined as "Homo sapiens" and the confidence value set as ≥ 0.4, followed by the PPI network analysis of the overlapping targets. e acquired data were imported into the Cytoscape software to analyze the network topology parameter value, thus screening the top ten targets of the degree value (DC value) as the key targets in the network.

Pathway Enrichment Analysis.
e overlapping targets obtained by the PPI network analysis were uploaded to the database for annotation, visualization, and integrated discovery (David) [20], followed by gene ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses with the condition value as p < 0.05. e top 10 pathways in the enrichment analysis were selected to draw a bubble chart.

Verification of Molecular Docking between Key Targets and Pharmaceutical Ingredients.
e key targets (top 10 target proteins) were selected in accordance with the median value of the PPI network.
ereafter, the corresponding chemical active ingredients related to these target proteins were identified in the protein interaction network to confirm the binding activity of these key target proteins to the chemical active ingredients. e TCMSP database was applied to acquire the chemical structure of the chemical active ingredients. e structure files of the target protein were identified in the RCSB PDB database (http://www.Rcsb.org/) with the species restricted to "Homo sapiens." We preferred the proteins with original ligand small molecules and high similarity between small molecular structure and the active ingredients of the docking chemical.
e Autodock Vina was employed to perform modification operations, such as dehydration and hydrogenation, on the retrieved target protein. e protein was output in PDBQT format, and molecular docking was implemented to determine the global optimal binding conformation according to the value of docking binding energy. PyMOL was utilized to visualize the docking result.

Statistical Analysis.
e data were analyzed using the GraphPad Prism 8 software. Continuous variables were analyzed by t-test. p < 0.05 was statistically different.

Changes in Immune-Inflammatory Indexes after XFC Treatment.
To observe the situation of the patients with AS after XFC treatment, we selected clinical immune-inflammatory indexes for validation, including ESR, CRP, IgA, IgM, IgG, C3, and C4. Compared with the indexes before treatment, the levels of ESR, CRP, IgA, IgG, C3, and C4 were significantly decreased (p < 0.01; Table 1), while there was no significant change in IgM.

Analysis of Association Rules between XFC and Immune-Inflammatory Indexes.
e minimum confidence degree was set to 50%, and the minimum support degree was set to 10%. e correlation between XFC and immune-inflammatory indexes was acquired using the Apriori module. Association rule analysis showed that the confidence degree of XFC with improvement on CRP was 78.23%. e confidence degree with improvement on ESR was 63.79%, and the confidence degree with improvement on IgA was 66.67%. Meanwhile, all the lifting degrees are greater than 1. ey indicated that XFC was strongly associated with the improvement of AS immune-inflammatory indexes. e results were shown in Table 2.

Potential Active Ingredients of XFC Drugs.
A total of 92 chemical active ingredients in the XFC formula were harvested from the TCMSP database, among which 20 were Astragalus, 9 were coix seed, 51 were Tripterygium wilfordii, and 12 were centipede (Supplementary Table S1).

e Prediction of the Targets of XFC-Active Ingredients.
e target prediction function of the TCMSP database was utilized. Subsequent to data screening and duplicate value removal, 208 predicted targets of chemical active ingredients in the XFC formula were yielded. e Cytoscape software was employed to complete the relationship network diagram of the effective chemical ingredients and predict the targets of the XFC formula ( Figure 1). Among them, there were 818 interaction relationships and 263 targets for active ingredients and targets.

e Disease-Related Targets of AS.
e targets obtained in each database were intersected, screened, and processed through the GeneCards database, which yielded 629 disease-related targets of AS.
e Venn diagram of "drugdisease-related target" (Figure 2) was plotted, which acquired 57 targets of XFC for AS treatment (Supplementary  Table S2).

PPI Network Analysis Results.
e PPI network diagram of XFC in AS treatment ( Figure 3) had 56 nodes and 704 interaction relationships. e target proteins were ranked based on the degree value (DC value), and the top 10 target proteins were selected, including IL-6, TNF, VEGFA, CXCL8, IL-1B, PTGS2, STAT3, CCL2, IL-10, and IL-4. ese 10 proteins were defined as the key targets in the PPI network diagram (Table 3).

Pathway Enrichment Analysis Results.
e enrichment analysis of overlapping targets (57 targets of XFC for AS treatment) was implemented using the David database, demonstrating that the overall treatment process had multiple signaling pathways that participated in the orchestration of various biological processes. e results of the GO function enrichment analysis displayed that the overlapping targets of XFC in AS treatment were mainly enriched in biological processes, such as inflammatory response, cell response to lipopolysaccharide, positive regulation of NF-κB transcription factor activity, and immune response. ese overlapping targets were majorly enriched in cell components, including extracellular space, extracellular zone, outer plasma membrane, plasma membrane, cell surface, nuclear chromatin, membrane raft, extracellular  Evidence-Based Complementary and Alternative Medicine 3 exosomes, neuron projection, and protein extracellular matrix. Additionally, these overlapping targets were mostly enriched in molecular functions, such as cytokine activity, transcription activator activity, RNA polymerase II core promoter proximal region sequence-specific binding, and protein binding (Figure 4). KEGG enrichment analysis manifested that the related targets were significantly enriched in signaling pathways, such as TNF signaling pathway, NF-κB signaling pathway, and cytokine-cytokine receptor interaction ( Figure 5).

e Results of Molecular Docking.
e key targets in the PPI network were selected, and the relevant components in reverse were identified in the light of the protein interaction network for molecular docking. When the binding energy between the ligand and the receptor was lower, the binding conformation was more stable and the possibility of interaction was greater. e validation of the molecular docking results documented that the binding energy of the effective chemical active ingredients to the key target proteins was less than −5 kcal/mol (Figure 6), suggesting that the binding activity of the compound to the receptor was stable. e   Evidence-Based Complementary and Alternative Medicine PyMOL software was applied to visualize the protein docking pattern of the top four target proteins with binding energy (<-9 kcal/mol) and the anti-inflammatory active ingredients (Figure 7).

Discussion
ere was no record of the disease name in ancient medical records about AS, and AS belongs to the categories of "Bi syndrome", "Gubi", "Shenbi", and "Dalv" according to the clinical manifestations [21]. AS mainly results from the invasion of external pathogens, such as wind, cold, and dampness. At the same time, the insufficiency in the congenital endowment or the acquired loss of nutrition results in the deficiency of spleen Qi in patients, and thus the body system is imbalanced because of insufficient righteousness to resist external pathogens to weaken the function of spleen deficiency and dampness. AS contributes to the mutual accumulation of phlegm and blood stasis, unclear blood stasis, back pain, and spinal rigidity. Hence, TCM believes that the pathogenesis of AS is attributable to the deficiency of spleen Qi and blood stasis [22]. XFC has been documented to assume a good therapeutic role in the treatment of AS. However, because of the numerous chemical active ingredients in TCM and the lack of molecular research, the action mechanism of TCM in the treatment process remains unclear. Based on network pharmacology, bioinformatics analysis, and molecular docking technology, this study analyzed the active ingredients, the potential targets, and the pathways of XFC in the treatment of AS to study its molecular mechanism.

Evidence-Based Complementary and Alternative Medicine
Formononetin is an effective ingredient in the legume red clover, which can effectively restrain the inflammatory response in the pathological process of numerous diseases [29]. Inflammation is the key factor leading to the clinical symptoms of AS, such as joint damage, disability, and comorbidities. Hence, anti-inflammation is the primary treatment strategy [30]. e major pharmacological effects of the active ingredients are relatively concordant with the anti-inflammatory mechanism of AS. e results of the GO function enrichment analysis reflected that the treatment process was highly correlated with inflammation, cell response to lipopolysaccharide, positive regulation of NF-κB transcription factor activity, and immune response. Among them, the inflammatory response is strongly associated with the pathogenesis of AS.
In the results of the KEGG pathway enrichment analysis, the XFC formula exerted anti-inflammatory effects in AS by manipulating the TNF signaling pathway, NF-κB signaling pathway, and cytokine-cytokine receptor interaction. TNF in the TNF signaling pathway and NF-κB signaling pathway is one of the cytokines in the acute phase of inflammation, which is implicated in the systemic inflammatory response, and it accelerates the production of interleukins. TNF can repress chondrocyte apoptosis and cartilage bone destruction by inactivating TNF-binding receptors to reduce matrix metalloproteinase [31,32]. PTGS2 is rich in the TNF and NF-κB signaling pathways. e official website of KEGG unraveled that PTGS2 in these two signaling pathways modulated inflammation. PTGS2 upregulation can lead to the proliferation of fibroblast-like synovial cells, which influences the levels of TNF-α, IL-1β, and IL-6 in the serum and aggravates the inflammatory response [33]. Moreover, triptolide selectively suppresses COX-2 expression in human synovial fibroblasts, and the production of PGE2 inhibits IL-1α and diminishes TNF-α expression in the synovium of rats with adjuvant arthritis induced by collagen [34,35]. e results of molecular docking unveiled that the docking energy of the key targets in the PPI network and their corresponding chemical active ingredients were all less than −5 kcal/mol, indicating that the chemical active ingredients in the XFC formula and the disease-related targets of AS produced better protein-binding activity. In addition,  Evidence-Based Complementary and Alternative Medicine the binding energies of quercetin, kaempferol, and triptolide to TNF-α and of formononetin to PTGS2 molecule docking did not exceed -9 kcal/mol. is result illustrated that the chemical active ingredients of the XFC formulas had excellent binding activity to inflammatory targets.
is article still has certain limitations. is article still lacks experimental verification although the mechanism behind the action of Chinese herbal compounds in the treatment of AS was preliminarily predicted and molecular docking verification was carried out. Evidence-Based Complementary and Alternative Medicine

Conclusions
is research was conducted based on the methodologies of network pharmacology and molecular docking. rough the evaluation of the pharmacokinetic parameters of the XFC formula components and the enrichment analyses, it was noted that XFC could orchestrate cytokine activity, transcription activator activity, RNA polymerase II core promoter proximal region sequence-specific binding, and protein binding. e molecular mechanism of XFC in the treatment of AS was mainly achieved by the active ingredients, including quercetin, kaempferol, and triptolide. Formononetin reduced the expression of inflammatory factors (TNF and PTGS2) via TNF and NF-κB signaling pathways, thereby restricting the inflammatory response in AS.

Data Availability
All relevant data are included in the manuscript and supplementary material.

Conflicts of Interest
e authors declare that they have no conflicts of interest.

Authors' Contributions
FYY and LJ contributed to the study design. FYY contributed to data analysis and wrote the first draft. WJT, XL, and HD contributed to the specimen and data collection. LJ, GJC, and XL supervised the project and contributed to manuscript revision. All authors reviewed and accepted the content of the final manuscript.

Supplementary Materials
Active ingredients of XFC (Supplementary Table S1) and 57 targets of XFC for AS treatment (Supplementary Table S2) are in Supplemental Files. (Supplementary Materials)